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Abstract 

We construct an extension of the Standard Model (SM) which is based on grand 
unihcation with Pati-Salam symmetry. The setup is supplemented with the idea of 
spontaneous flavour symmetry breaking which is mediated through flavon fields with 
renormalizable couplings to new heavy fermions. While we argue that the new gauge 
bosons in this approach can be sufficiently heavy to be irrelevant at low energies, the 
fermionic partners of the SM quarks, in particular those for the third generation, 
can be relatively light and provide new sources of flavour violation. The size of the 
effects is constrained by the observed values of the SM Yukawa matrices, but in 
a way that is different from the standard minimal-flavour violation approach. We 
determine characteristic deviations from the SM that could eventually be observed 
in future precision measurements. 
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1 Introduction 


Grand Unified Theories (GUTs) represent an attractive direction to search for extensions 
of the Standard Model (SM) of particle physics (see e.g. [^). The pecnliar strnctnre of 
the SM gange gronp and its fermionic representations observed at low energies hnds a 
more natnral-looking explanation when embedded into larger (bnt simpler) gronps that 
may be realized at higher energies. In particnlar, each generation of SM qnarks and 
leptons hts into a 16-dimensional spinor representation of the gronp S'O(IO) which 
then also contains a right-chiral nentrino that can be nsed to explain the tiny bnt non¬ 
vanishing nentrino masses via the seesaw mechanism [3]-[^. Theoretical complication 
re-enters when it comes to the qnestion of how the GUT symmetry is broken down to 
the SU{3) X SU{2) x U{1) symmetry of the SM. Different Higgs representations of the 
GUT gronp can be chosen, and the symmetry breaking fnrther depends on the parameters 
controlling the gronnd state of the scalar potential. It tnrns ont that the SM gronp can 
be reached from S'(9(10) via basically two separate rontes: either via the Georgi-Glashow 
gronp SU{5) [^, or via the Pati-Salam (PS) gronp SU{4:) x SU{2) x SU{2y (^[^. 

Adopting a non-snpersymmetric framework, we will focns on the PS gange gronp with 
an explicit Z 2 symmetry relating the conplings and representations of the two SU (2) gronp 
factors (which affect the left- and right-chiral SM fermions, respectively). Recent work 
in this framework 10-^ snggests that it is worthwhile to consider PS scenarios with 
an extended held content (as compared to minimal setnps) where the breaking from the 
PS gronp to the SM gronp involves several distinct mass scales which eventnally realize 
gange-conpling nnihcation at the GUT scale. Here we investigate a particnlar constrnction 
which also allows ns to address qnestions relating to havonr and havonr mixing. 

The idea of combining grand nnihcation with havonr is itself not new (for reviews 
18 2^). The intimate connection of the SM fermions, which is dictated by 


see e.g. 

the GUT mnltiplet strnctnre, generally relates the Ynkawa matrices of diherent sectors 
in an nnrealistic way. In order to accommodate the observed Ynkawa strnctnres, it is 
therefore necessary to introdnce more than one GUT-breaking Higgs held 24,25 . Flavonr 


symmetries are then imposed to control the resnlting Ynkawa patterns by snitably chosen 
vacnnm expectation valnes (VEVs) of scalar havon helds. Typically, the so-constrncted 
Ynkawa matrices arise as sums of mostly non-renormalizable terms 
constructions compatible with the PS gauge group, see e.g. 


29 38 


26 28 . For explicit 


In contrast to this traditional approach, we implement an idea hrst proposed by Grin- 

where the concept of gauged havonr symmetries 


stein, Redi and Villadoro (GRV) in 39 


is realized within a renormalizable held theory. Among others, this requires to introduce 
new heavy partners for the SM fermions in order to cancel gauge anomalies associated 
with the chiral representations of the havonr symmetry group (see also [^). The new 
fermions can then be used to mediate havour-symmetry breaking — which is realized by 
VEVs of new matrix-valued scalar havon held^— to the SM sector. The original GRV 
paper focused on the SM quark multiplets. With a minimal set of new fermions and 
scalars, the authors have shown that the SM Yukawa matrices Yu,d for up- and down-type 


Mhe possibility of obtaining hierarchical structures for flavon VEVs (and thus realizing sequential 


discussed, for instance, in 42 -50 . 


flavour-symmetry breaking as described in 41 ) by an appropriate flavour-invariant potential has been 
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quarks are effectively given by the inverse of the corresponding flavon VEVs. Here, the 
partner of the top quark t' plays a special role: the large Yukawa coupling of the top 
quark implies that t' can be relatively light (i.e. within the reach of direct searches at the 
LHC at CERN), and that it can have a substantial mixing with the SM top quark. This 
leads to interesting new physics effects in electroweak precision observables and flavour 
transitions involving the third generation of quarks. At the same time, indirect flavour 
bounds from the hrst and second generation (most notably from K — K mixing and rare 
kaon decays) are naturally satished (see also (^). Variants of the GRV idea in the con¬ 
text of the SM 


52 , GUTs 53-55 , and supersymmetric theories 56,57 have also been 


discussed in the literature. 

The new essential ingredients in this work are a consequence of embedding the GRV 
mechanism in an explicit left-right (i.e. Z 2 ) symmetric PS GUT with a Higgs bi-doublet. 
First, this requires to introduce (at least) two partners for each conventional PS fermion 
multiplet. Second, in order to obtain realistic masses and mixing for quarks and leptons, 
one has to introduce different types of flavon helds. In the minimal setup that we will 
construct for the quark sector, we consider flavons that transform as singlets and triplets 
under the SU{2) factors of the PS symmetry group. The breaking of the flavour symmetry 
by the singlet VEV alone treats up- and down-type quarks in the same way, such that the 
corresponding Yukawa matrices are identical, Y„ = Y^. An additional VEV for the SU{2y 
triplet flavon then guarantees that 7 ^ Y^. We also discuss how this framework can be 
extended to yield realistic masses and mixings for the charged leptons and neutrinos. 

Our paper is organized as follows. In Section we introduce the general setup of 
our PS flavour theory and discuss the various scales of symmetry breaking. Moreover, 
we deduce approximate expressions for the Yukawa matrices as well as the effective light 
neutrino mass matrix. Section focuses on the flavour structure of the quark sector. We 
derive analytic formulas for the transformations from the flavour to the mass basis in order 
to extract the (non-standard) couplings of the quarks to the SM gauge bosons and the 
Higgs. In Section]^ we describe the method we have used to scan over the parameter space 
of the model. The results of this numerical scan are presented and discussed in detail in 
Section We conclude in Section Supplementary material and more technical details 
are provided in the Appendix. 


2 Pati-Salam model with gauged flavour symmetry 

As outlined in the introduction, our starting point is a non-supersymmetric high-energy 
theory with an underlying Pati-Salam gauge symmetry. A manifest left-right symmetric 
setup is realized by enforcing a discrete Z 2 symmetry. The flavour symmetry group in PS 
is dehned by the independent transformations of the two fermion representations which 
are used to embed the SM fermions. In the following, we restrict ourselves to the flavour 
group SU{3)i X SU{3)ii, where the left-chiral SM fermions transform as triplets under 
SU{3)i while the right-chiral SM fermions furnish triplet representations of SU{3)ii. This 
flavour symmetry, which we assume to be gauged, is spontaneously broken by the VEVs 
of flavon helds. By the usual Higgs mechanism, the associated Goldstone modes will 
become the longitudinal degrees of freedom of sixteen new heavy gauge bosons of the 
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flavour groupj^ Notice that the maximal flavour-symmetry group in PS is smaller than in 
the SM, where we encounter a U{3)^ symmetry in the quark sector (as discussed in the 
original GRV paper [^) and a f/(3)^ symmetry in the lepton sector. To summarize, the 
full symmetry of the Lagrangian is given by 

g = (^SU{4) X SU{2) X SU{2y^ X (^SU{3)i X SU{3)n) x Z 2 , (2.1) 

'-V-^ '-V-' 

Pati-Salam flavour 

where the global Z 2 symmetry is realized such that a general representation O of ^ is 
transformed into O of ^ as 

Q = {LJc,oJ,oj'){uJi,un) —G = (chc,a;',a;)(a;77-,cJ/) . (2.2) 


Here and in the following, u denotes a representation of the corresponding group factor 
and uJ is its complex conjugate. As we will only introduce (pseudo-)real representations 
of SU{2) X SU{2y, we have dropped the complex conjugation of the corresponding rep¬ 
resentations in n on the right-hand side of Eq. (2.2). 


2.1 Particle content and renormalizable Lagrangian 


The left- and right-chiral SM fermions and (including the right-chiral neutrino) are 
embedded in the ( 4 , 2 , 1 ) ( 3 , 1 ) and ( 4 , 1 , 2 ) ( 1 , 3 ) representations of Q. Notice that the trans¬ 
formation of Eq. (2.2) maps the right-chiral fermions qn into the complex conjugate of 
the left-chiral fermions. Therefore, the chirality of the spinors remains unchanged under a 
Z 2 transformation. The SM Higgs is embedded in the bi-doublet H ~ ( 1 , 2 , 2 )( 1 , 1 ), which 
effectively generates the structure of a two-Higgs-doublet model. In order to generalize 
the idea of GRV, it is necessary to introduce a set of new fermions, denoted as 

The former, i.e. and transform as ( 4 , 1 , 2 )( 1 , 3 ) and ( 4 , 1 , 2 )( 3 , 1 ), while the 
complex conjugate of the latter, i.e. correspond to the Z 2 conjugated representations 
of Since we impose an exact Z 2 symmetry, we have to introduce both pairs. 

Due to the single Higgs bi-doublet il, a flavour symmetry breaking flavon field S 
transforming trivially under PS cannot discriminate between the up-type and the down- 
type sector. In contrast to the original GRV model 39 , it is therefore impossible to 


generate a non-trivial flavour structure without mediating the S'f/(2)' breaking of the PS 
sector to the flavour sector. As a consequence, we are led to introduce PS non-singlet 
flavon helds, where the simplest case is realized by T ~ ( 1 , 3 , 1 )( 3 , 3 ) and T' ~ ( 1 , 1 , 3 )( 3 , 3 ), 
in addition to the singlet flavon S. 

The particle content defined up to this point treats quarks and leptons on equal foot¬ 
ing. In particular, neutrinos would be Dirac particles which copy the hierarchical pattern 
of the up-type quarks. In order to account for the evident difference of quark and lep¬ 
ton flavour, we take advantage of the possibility to generate Majorana mass terms for 


^The Yukawa sector is additionally invariant under three independent C/(l) symmetries: two associated 
with SU{3)i X SU{3)ii, while the third relates to the neutrino extension of the theory. All three are 
broken spontaneously by flavons of the neutrino sector. Massless Goldstone modes can, however, be 
avoided by breaking the U{1) symmetries explicitly in the scalar potential, for instance by terms involving 
the determinant of flavon fields. 
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Pati-Salam Symmetry 

Flavour Symmetry 

VFV 


^f/(4) X SU{2) X SU{2y 

SU{3)i X SU{3)ii 


Ql 

( 4 , 2 , 1 ) 

( 3 , 1 ) 

— 

Qr 

( 4 , 1 , 2 ) 

( 1 . 3 ) 

— 

H 

( 1 , 2 , 2 ) 

( 1 . 1 ) 

Vu,Vd 

Sz 

( 4 , 1 , 2 ) 

( 1 . 3 ) 

— 

St? 

( 4 , 1 , 2 ) 

( 3 , 1 ) 

— 


( 4 , 2 , 1 ) 

( 1 . 3 ) 

— 

—H 

( 4 , 2 , 1 ) 

( 3 , 1 ) 

— 

T 

( 1 , 3 , 1 ) 

( 3 , 3 ) 

0 

r 

( 1 , 1 , 3 ) 

( 3 , 3 ) 

±t' M 

s 

( 1 , 1 , 1 ) 

( 3 , 3 ) 

sM 

©L 

( 1 , 1 , 1 ) 

( 5 , 8 ) 

— 

Or 

( 1 , 1 , 1 ) 

( 8 , 3 ) 

— 

Su 

( 1 , 1 , 1 ) 

( 6 , 1 ) 

SuK 

SI 

( 1 , 1 , 1 ) 

( 1 . 6 ) 

S'yK 

$ 

( 4 , 2 , 1 ) 

( 8 . 1 ) 

0 


( 4 , 1 , 2 ) 

( 1 . 8 ) 



Table 1: The particle content of the theory with imposed Pati-Salam and flavour symme¬ 
try. Left- and right-chiral fermions 4^l,r ai'e denoted by subscripts L and R, respectively. 
The VEVs of the scalar helds (no subscript L or R) are given in the rightmost column. The 
lower part of the table shows fields necessary for generating Majorana neutrino masses. 


electrically neutral leptons. Demanding renormalizability, it is necessary to extend the 
particle content by further fermionic and s calar helds. Although not unique, there exists 
a preferred extension (see Appendix A.l) in which additional PS-neutral fermions ©z, 
and Qn acquire Majorana masses via their coupling to new havour symmetry breaking 
havons S'y and S'^. Furthermore, ©z and S/j (©^^ and El) are coupled to one another by 
means of yet another havon held <h' (<h). Transforming as ( 4 , 1 , 2 ) under PS, the vacuum 
expectation value of the havon breaks the Pati-Salam gauge symmetry in the direction 
of the SM singlet, thereby projecting out the SM neutral component of This in turn 
couples to the neutrino component of Ql (either directly or indirectly via Sz and qn) 
which eventually generates a Majorana mass for the left-chiral neutrino Ql- 

The particle content of the model required for the description of quark and lepton 
masses and mixing is summarized in Table together with the corresponding trans¬ 
formation properties under the underlying PS and havour symmetry. Note that we list 
all fermions using only right-chiral degrees of freedom, i.e. 'ipR and the Dirac adjoint ipL 
(instead of ■^z)- As each particle is accompanied by its Z 2 conjugate partner, the PS 
anomaly [S'f/(4)]^ vanishes by construction. Concerning the havour anomalies [S'?7(3)/]^ 
and [SU (3) 77 ]^, it is straightforward to check that the number of fermionic triplets matches 
the number of fermionic anti-triplets. Notice that we have chosen the octet representa- 
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tion in 02, and Qr in order to reflect the dimension of the PS representation of the other 
fermions. As the adjoint does not contribute to the SU{3) anomaly, both flavour anoma¬ 
lies vanish identically. Mixed anomalies do not exist due to the absence of U{1) factors 
in Q of Eq. (2.1). Hence the model presented in Table is free of any gauge anomaly. 


With the particle content specified in Table it is straightforward to formulate the 
most general renormalizable Yukawa Lagrangian which is invariant under the full symme¬ 
try Q. Collecting all terms which involve the heavy fermions Ql,r (required to generate 
light Majorana neutrino masses) in /lyuk) write 


C 


Yuk 


~ -^Yuk 


+ c 


Yuk ) 


(2.3) 


where 


-^Yuk ~ ^ QlH Si? -|- El {ks S + kt T') Sr + M T^lQr + h.c. 

+ A H qR-\- [ns S -|- kr T) ^r -\- M ^r h.c. , (2.4) 

and 


~ 0L Er +10R s, Ql + h.c. 

-|- .=,L $ 0R -|- 4 0^ Si 0R -|- h.c. 

+ 0L^t0R +h.c. . (2.5) 


In Eqs. (2.4) and (2.5), the terms of the first line are related to those of the second line 
through the action of the Z 2 symmetry. The term in the third line of Eq. (2.5) is its 
own Z 2 conjugate. Due to these relations, the Yukawa Lagrangian depends on only 
four independent parameters: three dimensionless couplings A, kr, kr, and one mass 
scale M. Redefining the phases of the fermions, it is clear that A and M can be chosen 
real and positive without loss of generality. As kr and kr appear only in combination 
with the flavon fields which will generally acquire complex VEVs, their phases can be kept 
arbitrary at this point. With regard to the neutrino Lagrangian Cyy^, we are at this point 
only interested in a qualitative analysis. Therefore, we have suppressed explicit coupling 
constants, which would again be related by the Z 2 symmetry, in Eq. (2.5). 


2.2 Symmetry breaking and mass scales 

The flavour structure of the SM fermions originates in the flavon fields. Assigning VEVs 
of order M or higher for the components of S and T', breaks the flavour symmetry 
SU[3)i X SU[3)ii above the electroweak scale v. The scalar held T, on the other hand, 
must not get a high scale VEV as it transforms as an SU{2) triplet. However, a possible 
mixing with the electroweak Higgs doublets would generically induce a VEV (T) below 
the electroweak scale at order v'^/M. As the impact of this VEV on the flavour structure 
is negligible compared to (S), we will approximate (T) ~ 0 in the following, for simplicity. 
In order to obtain a phenomenologically realistic difference between the effective Yukawa 
matrices of the up and the down quarks, we need the singular values of the VEVs (S) 
and (T') to be of similar size. 
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In addition to breaking the flavour symmetry, the VEV (T') breaks S'f/(2)' as well. In 
order to keep the standard charge assignments for up- and down-type quarks, we choose 
(T') to be aligned along the Ts direction of 517(2)'. Its smallest singular value is considered 
to be (at least) of the order of the scale M. Factoring out this explicit mass scale, we 
parameterize the vacuum structure responsible for the Yukawa matrices of the charged 
fermions as follows. 


Ks {S) = s M and kt {T') 


t' 0 
0 -t' 


M, 


(2.6) 


where s and t' are defined as dimensionless 3x3 matrices in flavour space. 

Turning to the neutrino sector, we adopt the double seesaw mechanism involving 
q^, Sr and ©r 0 The PS-neutral fermions 0 r acquire a Majorana mass (S,^) at a high 
scale Aty. Furthermore, a Dirac mass term 0 r(<I>')Sr is generated at the scale 
With this hierarchy of scales, the first stage of the double seesaw mechanism gives rise to 
a Majorana neutrino mass for the SM neutral component S^ within Sr. In the second 
stage of the double seesaw, light Majorana neutrino masses for the left-chiral neutrinos 
are induced via the Dirac coupling of q’l and S^. Analogous to Eq. (2.6), we parameterize 
the VEVs of the scalar fields required for generating light Majorana neutrino masses as 


{Su) = SuA^ 


{Sl) = slA, 


and (d)') = if' A^ 


(2.7) 


Here we have included the VEV of which will turn out to be crucial for the construction 
of the lepton sector as discussed in Appendix A.l The quantities Siy, s'^ and ip' are again 


dimensionless tensor structures in flavour space. Note that <1> transforms non-trivially 
under the electroweak symmetry. Similar to T, it must therefore not receive a high scale 
VEV, and we will approximate (<1>) ~ 0 in the following. 

In contrast to the original GRV model proposed in 


39 


in our setup the flavour gauge 
bosons do not play an essential role for the low-energy phenomenology. The reason for 
this is the inclusion of the lepton/neutrino sector which requires the introduction of the 
flavons S^, 5/ and d*' in addition to S and T'. While the latter acquire VEVs reaching 
down to the scale M (TeV regime), the former obtain their VEVs at much larger scales. 
Ay and A^. As these are related to the type I seesaw scale via A^/Aj, ~ Mgeesaw, it is 
reasonable to assume a typical hierarchy of scales. 


lO'^GeV < A^ < A, < Mpi^nck ~ GeV . 


18 , 


( 2 . 8 ) 


Here, the lower bound is relatively flexible (due to the vague definition of Mgeesaw) and 
has been chosen rather conservatively. As a consequence, the gauged flavour symmetry 
gets broken at a very large scale, and the associated flavour gauge bosons become too 
heavy to be relevant for current or future particle physics experiments. In principle, a 
subset of the flavour gauge bosons could remain massless down to lower scales. However, 
if the symmetric 3x3 matrix Sy has non-vanishing and non-degenerate singular values, 

^There exists another contribution to the effective light neutrino masses which additionally involves 
the fields qn, and Sr. Its structure will automatically become manifest when integrating out the 
heavy degrees of freedom in /lyuk- 
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10^ GeV IQio GeV GeV 


M S ,T' Mpgr ~ M^xm^ 

Figure 1: Illustration of the hierarchy of scales introduced in the PS model with gauged 
flavour symmetry. The numerical values are shown as an indicative example only. 
denotes the scale of the flavour gauge bosons, while Mpg- is the scale of the heavy PS 
gauge bosons. 


SU{3)j gets fully broken and all its flavour gauge bosons acquire masses at the scale A^. 
Similarly, all flavour gauge bosons of SU{3)ii become heavy due to (SI). 

Concerning the masses of the gauge bosons associated with the PS symmetry, we 
point out that we do not fully specify the breaking of S'17(4) or SU{2)'. Yet, already the 
flavon d)', which we have introduced for the purpose of generating light neutrino masses, 
breaks PS down to the SM 10,59 . Therefore, the masses of the non-SM gauge bosons, 
including the W'^ and the Z', are bounded from below by A^. Additional PS breaking 
helds may lead to even larger masses. 

As discussed above, the various scales of the model are only partially constrained. 
Their ordering is, however, determined as illustrated in Fig. In addition to the hierar¬ 
chies of scales, the singular values of (S') and (T') feature internal hierarchies, with the 
smallest reaching down to the scale M. The value of M is bounded from below by flavour 
precision observables as well as the non-observation of new charged fermionic states. For 
the purpose of our numerical parameter scan, we will assume values of a few TeV. Due 
to the strong hierarchy of quark masses, the largest singular values of (S') and (T') take 
numerical values up to order 10^° GeV. This is in fact the reason for choosing the lower 
bound on A^ as given in Eq. (2.8). On the other hand, A^ is bounded from above by the 
GUT scale Mqut, as (d)') breaks PS down to the SM. Finally, the largest scale Ay must 
lie between A^ and the Planck scale Mpianck- In addition to these scales. Fig. also shows 
the mass scale of the flavour gauge bosons as well as the mass scale Mps" of the 

heavy PS gauge bosons as discussed above. 

Strictly speaking, the model is formulated at the GUT scale. Due to the presence 
of multiple scales, it is generally necessary to consider threshold corrections at each of 
these scales. In what follows, we will, however, neglect such effects as we do not aim 
for a comprehensive top-down description of the model. We will therefore perform our 
numerical scan using flavour parameters (e.g. quark masses) which are hxed at the Z boson 
mass scale. 


2.3 Approximate flavour structure of quarks and leptons 

In a hrst step to understand the structure of the low-energy effective Lagrangian derived 
from Eq. (2.3), we formally integrate out the heavy degrees of freedom. The result will 
provide reasonably simple expressions for the approximate flavour structure of quarks 
and leptons. We begin our discussion with the quark sector, where all relevant terms are 
included in of Eq- (2-4). Assuming for the moment all VEVs of S and T' to be 
signihcantly larger than the intrinsic mass scale M, i.e. s ~ F S> 1, we can integrate out 
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(2.9) 


the heavy fermions and Dehning the 3x3 matrices 

= s +1' , td = s-t' , 
the resulting effective Yukawa potential takes the form 


^Yuk 


= ->^hlql [t^^ + s ql +h.c. 
-^h^dQi [td^ + 5 "^] 4 + h.c. , 


( 2 . 10 ) 


where the 517(2) and 517(2)' doublets have been explicitly decomposed into their up and 
down components. The effective couplings correspond to the standard quark Yukawa 
matrices in a two-Higgs-doublet model, and can be directly read off as 


^ -A + s ’] 


( 2 . 11 ) 


It is worth mentioning that this relation between the flavon VEVs and the SM Yukawa 
matrices is neither a linear one as in Minimal-Flavour Violation (MFV) [^-64 
simple inverse one as in 


39 


nor a 


where the difference with the latter can be traced back to 
the imposed Z 2 symmetry. Yet, similarly to the original GRV case, the expressions in 


Fq. (2.11) show that the observed hierarchy of the Yukawa matrices is turned into an 
inverse hierarchy of the 3x3 matrices s and t'. Fq. (2.11) can now be used to £t the 


19 independent parameters of s and t' (six singular values, six mixing angles and seven 
phases) to the quark masses and the CKM matrix Vckm- In contrast to the SM, the 
mixing of the right-chiral quarks (parameterized by in the following) is physical in 

our PS model and has to be kept as a consequence. As the relation between Yu^ and s,t' 
is given by a coupled matrix equation, it is a non-trivial task to invert Eq. (2.11). We 
describe the method used to obtain numerical results for s and t' in Appendix |B[ 

An important feature, already arising at this stage, is the existence of multiple so¬ 
lutions when inverting Eq. (2.11). This can be easily understood in the case of one 


generation, where the resulting two solutions can be expressed analytically. Explicitly, we 
End 


-8A 


s = 


tu 


^Vu + ^Vd T S-y /{Uu Vd) T 9 VuUd 

-8A 

^Uu ‘^Vd i 3 {qju Vd) T 9 VuUd 


(2.12a) 

(2.12b) 


with denoting the one-generation up-type and down-type Yukawa couplings. Using 
these expressions it is straightforward to determine t' as well as td- Extending this result 
to three generations gives rise to a total of eight solutions for s and t'. We emphasize that 
all of these generate the same approximate Yukawa matrices via Eq. (2.11). However, 


as we will see in Section different solutions may generate different phenomenological 
predictions beyond the SM. 


^To be precise, considering the quark sector, we only have to integrate out the colour (anti-)triplets 
contained in the (anti-)fundamental representations of SU^A). 
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Turning to the lepton sector, we first remark that the structure of the charged lepton 
Yukawa matrix is identical to the one of the down-type quarks. This is a common feature 
to the simplest Pati-Salam models and a reasonably good hrst approximation. Possible 
modihcations to the charged lepton sector are briefly sketched in Appendix |A.2[ In the 
neutrino sector, the Lagrangian involving the heavy fermions Ql,r generates a Majorana 
mass term for the neutrino component within the PS multiplet Integrating out 


in /lyuk of Eq. (2.5) yields 


C 


Yuk 


2 A,. 






(2.13) 


We remark that Sy is a symmetric 3x3 matrix in SU{3)i flavour space, while (p^A,^ 
denotes the VEV of <!>' which points in the direction of the SM neutral component of 
the PS representation ( 4 , 1 , 2 ). Being an octet of SU{3)jj, with a = 1,...8 does not 
affect the flavour structure of the neutrino sector. Combining this Majorana mass term 
for with the couplings of the electromagnetically neutral fermions in /lyuk of Eq. (2.4) 
generates Majorana masses for the left-chiral neutrinos. The resulting structure of the 
neutrino mass matrix can be obtained by formally integrating out the heavy neutral 
fermions and as well as and It is given by 


m 


eff 


A. 


2«A, 


(l + s s„ {l + s 


(2.14) 


where Vu/\/^ denotes the VEV of the neutral Higgs cf. Appendix A.l, It is important 


to note that the flavour structure of the light neutrinos is decoupled from the quark sector. 
In particular, the distinct hierarchical structure of s and tu is approximately cancelled in 
the combination s~^tu, such that the structure of the effective neutrino mass matrix is 
mainly determined by As a consequence, the PMNS mixing of the lepton sector can be 
signihcantly different from the CKM mixing of the quarks. Assuming the dimensionless 
parameters s~Hu, Sy and to be of order one, we hnd that Ay ~ 10^®GeV ~ Mgut 
and Ky ~ lO^^GeV Afpianck gives rise to realistic neutrino masses of the order 0.1 eV. 
As such, the model provides sufficient freedom to construct a phenomenologically viable 
lepton sector. More detailed studies are left for future work. 


3 Quark flavour sector 


The effective quark flavour structure derived in Section 2.3 is only an approximation which 
relies on the simplifying assumption that s,t' 3> 1. Glearly, such a relation is not satished 
for the Yukawa couplings of the third generation of quarks which should be of order one 
for the top quark (and, depending on tan/3 = Vu/vd-, also for the bottom quark). Taking 
this fact into account, it is necessary to diagonalize the full 9x9 mass matrices of the 
SM quarks and their heavy partners. The corresponding transformations are then applied 
to the gauge-kinetic sector in order to extract the characteristic features of the effective 
low-energy quark flavour sector. We also comment on the coupling of the Higgs held to 
the SM quarks, which is no longer fully aligned with the diagonal quark mass matrices. 
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3.1 Diagonalizing the quark mass matrices 


In this subsection, we re-examine the Yukawa Lagrangian of Eq. (2.4). Assigning VEVs to 


the flavour symmetry breaking helds S and T' as well as the SU (2) x SU (2)' bi-doublet H, 


Eq. (2.4) gives rise to 


A + M + M S; + h.c. 

+ ^ A M s M h.c. 

A (rt chj . 


(3.1a) 

(3.1b) 


Collecting the left-chiral (right-chiral) fermions in 4//^ these bilinear mass terms 

can be compactly written as and For definiteness, we focus on the 

up-type quark sector in the rest of this subsection. Analogous results arise for the down- 
type quark sector, and the corresponding expressions are obtained by simply replacing 
the index u by d. Defining 


A^K,E“,Hp, ^ (5S,E“,Ha 


(3,2) 


the 9x9 mass matrix takes the form 


= 



(3.3) 


where each entry corresponds itself to a 3 x 3 matrix in generation space. The quantity 


^ 7 /. - 


V2M 


< 1 


(3.4) 


can be treated as a small expansion parameter. In the following, we describe the sequence 
of basis transformations which diagonalizes of Eq. (3.3) in powers of up to quadratic 


order. More technical details of this diagonalization are relegated to Appendix 

In the hrst step, we make use of the SU (3)/ x SU (3)// flavour symmetry to diagonalize 
the 3x3 matrix s in Eq. (3.3). This corresponds to a choice of basis which can always 


be made without loss of generality. We can thus simply replace s by s, where here and 
in the following the hat denotes a diagonal matrix. Next, we apply a bi-unitary basis 
transformation to diagonalize the submatrix tu, 


i.e. 


L - V,, t„ u\. 


(3.6) 


Having exhausted the flavour symmetry by choosing a diagonal s, the unitary 3x3 
matrices W and Uu will reappear elsewhere in the full mass matrix. As discussed in 
Appendix they can however be shifted to the e„-suppressed blocks. In the limiting case 
of e„ = 0, the resulting mass matrix takes the form of Eq. (3.3) with s and tu replaced 
by s and tu. In this limit, the three generations decouple from one another and the 
9x9 matrix decomposes into three 3x3 blocks, one for each generation i. Introducing 
only two mixing angles for each generation (whose values depend on s and tu), it is 
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straightforward to diagonalize these blocks exactly. In the third step, we therefore apply 
such a transformation to the full mass matrix including the e„-suppressed blocks. This 
gives rise to a mass matrix of the form 





(3.6) 


where the 3x3 matrices a„, bu, Cu, du e„ and / depend on I 4 , as well as the above 
mentioned mixing. The explicit expressions are given in Appendix 0 Note that the 
mass matrix in Eq. (3.6) is diagonal at zeroth order in the expansion parameter. Block- 


diagonalization to second order in requires an intricate basis transformation which can 
be found in Appendix Yet, the resulting mass matrix simply reads 





0 
0 

f + 0{el), 


M + 0{el)M 


( 3 . 7 ) 


The remaining transformation needed to render this matrix completely diagonal (up to 
order e^) involves only rotations within the three non-vanishing 3x3 blocks. These do 
not mix the light quarks with the heavy partners. Due to our lack of experimental data 
on any heavy fermions, it is safe to ignore the corresponding 3x3 rotations which would 
take the general form 1 -|- 0{e^. For the three light quarks on the other hand, we dehne 
the bi-unitary transformation which diagonalizes a^. 


K = V,,a„Ul 


(3,8) 


The complete sequence of transformations diagonalizes the 9x9 quark mass matrix of 
Eq. (3.3) up to second order in e^. It therefore constitutes the change from the original 
flavour basis to the approximate mass basis The individual steps of this change 

of basis are explicitly given in Appendix In the next subsection, we will apply these 
transformations to the gauge-kinetic sector. 


3.2 Gauge-kinetic terms 


As discussed in Section |2.2[ the flavour gauge bosons are far beyond the reach of current 
or future particle physics experiments. They are therefore phenomenologically irrelevant 
and we do not consider them any further. Concerning the 15 gauge bosons associated 
with the S'17(4) factor of the PS symmetry we note that 6 become extremely massive and 
irrelevant when SU (4) —)■ SU ( 3 )c x f/(l) b-l- Moreover, the 8 gluons of SU ( 3 )c are flavour 
blind and can be ignored as well. The remaining U{1)b-l boson will mix with the neutral 
gauge bosons of SU (2) x SU (2)' and has to be kept. Being mainly interested in the flavour 
structure, it is therefore sufficient to focus on the gauge-kinetic terms corresponding to 
the SU{2) X SU{2y x U{1)b-l part of the PS symmetry. In order to be able to apply 
the sequence of basis transformations described above (and in Appendix]^, we express 
the relevant gauge-kinetic terms by means of the left- and right-chiral vectors ^i and T/j 
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defined in Eq. (3.2). In this context, it is important to realize that does not solely 
contain donblets of SU{2) bnt also donblets of 517(2)'. Likewise, d'R contains donblets of 
both SU{2y and 517(2). Taking this fact into acconnt, the relevant gange-kinetic terms 
read 


Tikin D 

+ (^7#'r) 

+ ^l{19b-lQb-l$b-l^'^l + 9b-lQb-l$b-l^ ^ (3-9) 

where Qb-l is the difference of baryon and lepton nnmber, and the /C matrices 


/I 0 0\ 
/Cl = 0 0 0 

\0 0 1 / 


/O 0 0\ 

/Ck = 0 0 0 , (3.10a) 

\0 0 1 / 


/O 0 0\ 
/Ci = 0 1 0 
\0 0 0 / 


/l 0 0\ 

/c^ = 0 1 0 

\0 0 0 / 


(3.10b) 


encode the non-trivial 517(2) x 517(2)' transformation properties of Tl and Here we 
emphasize that the 517(2) gange bosons conple to both left- and right-chiral fermions. 
Likewise, the 517(2)' gange bosons conple to fermions of both chiralities. 

While breaking 517(2)' x 17(1)l_l —t 17(l)y the nentral bosons and Bb-l 

mix, resnlting in a massless By and a massive Z' boson. Analogons to the electroweak 
symmetry breaking in the SM, this mixing is described by s'^, and c'^ which satisfy 
9 s'w = Qb-lc'-^y = gy. The charged 517(2)' gange bosons W'^ = =F be¬ 

come massive from the 517(2)' breaking as well. As this breaking arises at an extremely 
high scale, cf. Fignre[^ the W'^ and Z' gange bosons are irrelevant for low-energy flavonr 
effects. In a second step, the breaking 517(2) x 17(l)y — )■ 17(l)em at the electroweak scale 

indnces a mixing of the nentral gange bosons and By which is parameterized by sw 

and cw satisfying g sw = gv cw = e. Using these relations and the two-step mixing 

$B-L = sw$) - , (3.11a) 

= Sy^{cwA- sw$) + , (3.11b) 

= Sw/A + c-wA 1 (3.11c) 


we can rewrite the gange-kinetic Lagrangian, keeping only the terms which are relevant 
for low-energy flavonr effects, 

/ikin 3 

+ ({(AvK^L — S^^'l)t^ — \s^Qb-i\A^L 

Cw \ J 

+ ^Le{{lCL + lCy)T^ + \QB-i^A^L 

+ {L^ R) . (3.12) 
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In order to simplify the notation, it is convenient to decompose the isospin doublets in 
the charged current interactions explicitly into their components. Using the relations 
K^l + = 1 as well a^ Qe= \ Qb-l + we obtain 

Ain 3 -^Wl>CLM'*1>i + h.C. 

+ — ((t" - o'ivQ,) 1 - ic'l a ^*1 

Cw \ J 

+ — -si, Qel + }Cr tA 

Cw \ J 

+ e(^^Qe4^L+ "^nQeA^R) . (3.13) 


When applying the sequence of basis transformations discussed above, it is important to 
take into account the isospin structure of the charged and neutral couplings. As a result, 
the flavour structure of all terms involving the /C matrices is modihed in a non-trivial 
way, while the form of the remaining terms is left unchanged. For instance, the last line 
of Eq. (3.13) shows that the coupling of the photon is always diagonal and proportional 
to the corresponding electric charge. The explicit expressions of the full /C matrices in 
the approximate mass basis are provided in Appendix They will be important when 
studying the gauge interactions of the light quarks with their heavy partners such as 
e.g. the coupling bt' W~. Being mainly interested in the phenomenological flavour effects 
involving only the three generations of light quarks ^ (where the prime denotes the 
approximate mass basis), we focus our attention on the upper left 3x3 block. Separating 
the couplings which are already present in the SM (such as Vckm) from terms which are 
characteristic to our setup (i.e. AVckm, Uckm, ^Qzqj^q^ and Agzq^q^^ fhe terms of the 
gauge-kinetic Lagrangian take the form 


■^kin 


+ 


q'l (^ckm — AVckm) + A-c. 

Ar Uckm Ar + h.c. 

^ - SwQe) 1 - ^ 9 zq^q^ t Ql 

-SwQel + AgZq^q^A'^^q'q^ . 


(3.14) 


We emphasize that the two matrices Agzq^^q^ and Agzq^^q^^ will generally be different for 
the two isospin components. Mindful of this subtlety, we however suppress the corre¬ 
sponding isospin indices for the sake of notational simplicity. Comparing the explicit 
expressions of Appendix [C| in particular Eq. (C.ll), with Eq. (3.14) reveals that the 


®The relation between the electric charge Qe and Qb-l takes this simple form as the fermions of the 
model transform solely under one of the two SU{2) gauge factors. Therefore, is understood to act on 
either SU{2) or 811(2)', depending on the specific particle. 
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Figure 2: Leading order corrections to the quark masses and the Higgs-quark-quark cou¬ 
pling ghgq. 


low-energy flavour effects are parameterized as follows, 

hcKM = 1 (3.15a) 

AUkm = \ K ( 6 . bl 4 + b, ef 6 + , 5 ) yt + q , (3.15b) 

UcKM=Uuc\f~‘^ CdUd^ eued +O , (3.15c) 

/^gz-g^g^=Vhe-H^V^e^ + 0{e^) , (3.15d) 

= u cWt + O (e^) . (3.15e) 


We point out that the effective CKM matrix Id^KM is no longer unitary as it receives a 
correction AVckm- Like all the other beyond the SM flavour effects, this correction is of 
second order in the expansion parameters eu,d- 


3.3 Non-standard Higgs couplings 


As can be seen from the Lagrangian of Eq. (2.4), there exists no direct coupling of the 
Higgs boson to a pair of SM quarks. Similar to the quark masses, such an interaction 
is generated effectively through couplings to the heavy quark partners. For this reason, 
the Higgs-quark-quark coupling g^qq is no longer exactly diagonal in generation space nor 
proportional to the quark masses. The effective expression for ghqq can be deduced by 
formally integrating out the heavy fermions or, alternatively, by explicitly applying the 
sequence of basis transformations defined above to the Higgs couplings. As a result, we 
And that the corrections to the quark masses as well as the Higgs-quark-quark coupling 
are both of order e\ however, they differ from one another by a combinatorial factor of 
three. 

This effect can be understood by explicitly studying the leading order corrections 


starting from the basis of Eq. (3.6). We stress that this basis does not make use of any 


approximation in terms of a truncated eu,d expansion. Reinserting the original Higgs held 
in Eq. (3.6), the corresponding diagrams, depicted in Figure]^ lead to an effective coupling 
of three Higgs bosons to a pair of SM quarks. They provide a correction to the quark 
masses if all Higgs fields acquire their VEV. Alternatively, if only two of them are set to 
their VEV, the diagrams of Figure contribute to the Higgs-quark-quark coupling ghqq. 
With three possible ways of choosing two Higgs VEVs (and one Higgs held), the latter 
correction picks up the combinatorial factor of three. Such an effect also occurs in the 


SM when considering general effective dimension six operators 65 . 








As the effect arise identically in the up-type as well as the down-type quark sector, we 
suppress the subscripts u and d in our resulting expression which are valid up to order 
In the approximate mass basis, where the three light generations of quarks are denoted 
by q'j^ j:j, the effective Yukawa coupling takes the form 


Y ^ Y - 

2 L 


while the Higgs-quark-quark coupling is given by 

3e2 ^ 


Vbe-‘^b^V^ - Y + Y ■ U f-^ cU^ 


9hqq 


Y - 


Vbr'‘h‘<V' - Y + Y ■ Uc^ / ^ cU' 


(3.16) 


(3.17) 


In these expressions, Y is the exactly diagonalized version of the matrix a, using the 
unitary transformations V and W, cf. Eq. (3.8). We observe that the Yukawa matrix in 
Eq. (3.16) is only diagonal up to hrst order in eu^d- Diagonalizing Y to second (or in 
fact arbitrary) order can be achieved by small modihcations of V and U. Applying the 
same basis transformation to Qhqq of Eq. (3.17) will generally leave the Higgs-quark-quark 
coupling non-diagonal at order Moreover, comparing Eqs. (3.16) and (3.17), is it 
possible to show that the diagonal couplings of ghqq are reduced compared to their SM 
values. Numerically, we explicitly End such an effect in our scan over model parameters. 


where we exemplarily show the (normalized) Higgs-top-top coupling in Figure 7(a) 


4 Scanning the parameter space 


In the previous section, we have calculated analytical expressions for the Yukawa matrices 
and deduced flavour parameters of the model in terms of the underlying VEVs (S') and 
(T'). However, we are not able to solve the derived relations, such as e.g. Eq. (3.8), 


analytically for s and t'. In other words, we cannot parameterize the flavour signatures 
of the high-energy model using the physically known low-energy flavour parameters, i.e. 
the quark masses and the CKM mixing Vckm- 


On the other hand, it is possible to invert the approximate formula of Eq. (2.11) numer¬ 


ically. Doing so, we encounter the technical complication of multiple solutions, as already 
mentioned at the end of Section |2.3[ Moreover, starting from a left-right symmetric the¬ 


ory the right-chiral up-type and down-type quark sector cannot be rotated independently. 
This enforces the introduction of the mixing matrix D^km which parameterizes the cou¬ 
pling of the W'^ bosons to the right-chiral quarks and contains six physical phases. As 
both the mixing matrix U^^iM appropriate choice of a solution are unknown, 

we choose to scan over these degrees of freedom. In this section, we describe the technical 
details of our procedure. In principle a scan of the parameter space should be performed 
over all parameters of the model, namely s, t'. A, tan/3 = Vu/vd and M. However, as s 
and t' feature large hierarchies a direct scan over these is not reasonable. Thus we use the 
alternative ansatz and scan over so-called “adapted flavour parameters” which we define 
in the following. We point out that these are deduced from the parameters of the SM 
Yukawa couplings. 
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rhu 

0.5 - 

2.9 MeV 

Fid 

1.2 - 

4.8 MeV 

^12 

12.89° - 

13.19^ 

fhc 

0.53- 

0.71 GeV 

nis 

30- 

78 MeV 

^23 

1.54° - 

2.56° 

fht 

162 - 

288 GeV 

rhb 

2.78- 

4.44 GeV 

^13 

0.101° - 

- 0.28° 


Table 2: Ranges of the adapted flavour parameters over which the scan is performed. 
The mixing angles of Uqkm varied in their full range [0°,90°] and additionally in 
a reduced range [0°,1.5°]. All phases are varied between 0 and 27r. 


4.1 Adapted flavour parameters 


We start our considerations with the quark masses and mixings as htted in the SM 
framework [^. In order to take into account the physical but experimentally unknown 
mixing of the right-chiral quarks in the S't/(2)' charged current, we also add mixing angles 
and phases for the matrix ?7 ckm- With this input we calculate an explicit numerical 
realization for the 3x3 matrices W and in an arbitrary basis. For a given such pair it 
is possible to invert Eq. (2.11), thereby deriving s and t' numerically, see Appendix]^ As 
explained in Section 2^, these are however only approximations, and we expect sizeable 
corrections for the third generations. Nevertheless they provide a well-motivated starting 
point for the exact computation. 

Inserting the so-derived numerical values of s and t' into the 9x9 mass matrix of 
Eq. (3.3), we can diagonalize (and likewise explicitly. As our experimental 


knowledge is limited to the three generations of light quarks, it is sufficient to focus on 
the upper left 3x3 blocks of both the full mass matrices as well as the associated mixing 
matrices. Comparing these with the input parameters, we hnd that the masses and mixing 
angles related to the third generation are systematically to small. We therefore adjust 
the input parameters such that the effective 3x3 Yukawa matrices of the full theory 
match the matrices Yu and Yd of the SM. We denote these adjusted input parameters as 
“adapted flavour parameters” (labelled by a tilde), and dehne ranges (given in Table 
over which we perform the scan. 


4.2 Details of the scan 


Having dehned the adapted flavour parameters we describe in detail our systematic scan 
over possible flavour effects. We do not claim this scan to be exhaustive as the procedure 
may systematically exclude allowed regions of the underlying parameter space. Neverthe¬ 
less, it provides important insights into the essential flavour effects of this model. The 
scan itself is performed using the following steps. 


1. Randomly generate a point in the space of adapted flavour parameters where the 
allowed ranges are given in Table Furthermore, choose also random values for 
A e [1.5,3], tan/? e [1,15] and M e [750,2500] GeV. 


2. Calculate s and t' by inverting Eq. (2.11) using the procedure described in Ap¬ 
pendix Here we randomly choose one of the eight solutions. 
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3. Insert the so-derived s and t' as well as A, tan/3 and M into the fnll 9x9 mass 
matrices and diagonalize these nnmerically. 

4. Dednce the effective SM-like Yukawa parameters for the three light generations, i.e. 

and 

5. Compare these with the experimentally allowed 3cr ranges]^ In case of agreement, 
save the point (s, t'. A, tan /3, M) as a viable choice of input parameters. 


Following these steps, we have generated about 3,000 viable points using the full range of 
adapted parameters. However, scanning the full parameter space of the mixing angles 6'^ 
of is rather inefficient as only a small fraction of all points turns out to be physically 

acceptable. For this reason, we have additionally performed a more extensive scan where 
we choose random values of within the small interval [0°,1.5°]. The resulting 30,000 
viable points differ qualitatively from the full scan only in Uckm, i-e. the coupling of 
the right-chiral quarks to the bosons. In order to enhance the statistics, we have 
combined both sets of points in our analysis (unless we are interested in the results for 
h^CKM itself). 

Before turning to the presentation and discussion of the physical results of our scan, 
we comment on the qualitative structure of the flavon VEVs (S') and (T'). In the original 
GRV model 39 , the flavon VEVs can be derived from the SM Yukawa matrices via 
the approximate relation (S') cx: Y~^. This simple relation entails that the flavon VEVs 
directly inherit the hierarchies of the Yukawa matrices in an inverse way. In contrast to 
this, the situation in our PS flavour model is more involved. Due to the linear combination 
of the VEVs (S') and (T') in Eq. (2.11), a simple direct relation of the hierarchies in the 
SM Yukawas and the flavon VEVs does not exist. While many viable points of the scan 
do indeed feature flavon VEVs with hierarchies similar to those of the inverse Yukawas, 
there exists a signihcant number of points where the hierarchies of VEVs is less (or more) 
pronounced. Similarly, the mixing angles describing the mismatch of the matrices s and t' 
can in some cases be signihcantly larger than those of the CKM matrix. 


5 Quark flavour phenomenology 


In this section, we discuss the physical implications of our PS flavour model. Starting 
from the adapted flavour parameters together with the ranges for A, tan/3 and M as 
dehned in Section we observe that the values of tan /3 and M corresponding to the 
viable points are equally distributed throughout the full allowed ranges. In contrast to 
this, the parameter A does not cover the full allowed range but only A G [1.5,2.5]. In fact, 
it peaks around A = 2, which is exactly the value we have used to determine the ranges 
of the adapted flavour parameters. In principle, we could have expanded our scan by 
additionally deriving adapted flavour parameters for other values of A, thereby enlarging 


®We consider the masses at the electroweak scale 


• u.d 

I.e. m- 


{Mz) 


67 


and the experimental constraints 


on the absolute values of the individual CKM entries (not those obtained from additionally demanding 
unitarity of the CKM matrix). 
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the intervals in Table However, we refrain from doing so as we are mainly interested 
in this part of the parameter spacej^ 

Regarding the M-dependence of all physical resnlts, we note that this has been fac¬ 
tored ont in the dehnition of the flavonr breaking VEVs of Eq. (2.6). Following this 
parameterization, the approximate Yukawa relations of Eq. (2.11) do not explicitly de¬ 
pend on M. Likewise, we do not expect any explicit M-dependence of the 3x3 matrices 
a... f, cf. Eq. (3.6). On the other hand, an implicit M-dependence may in principle be 
induced via the intricate procedure of determining s and f', in particular for the third 
generation of quarks. We have therefore checked for such a possibility in our scan. How¬ 
ever, no extra (implicit) dependence on M is found in any of the quantities, at least above 
a threshold of M ~ 1 TeV. This observation allows us to understand the dependence on 
the scale M explicitly throughout the model, particularly in all corrections. 

As an application of this fact, we may limit our discussion of characteristic flavour 
effects to a narrower mass band. In some of the following studies we choose this to 
be M G [1,1.2] TeV, which reduces our set to approximately 3,500 viable points. This 
restriction clarihes the dependence of physical observables on the other input parameters 
of the scan as it reduces the spread due to the variation of M. The results obtained for 
this reduced mass band may afterwards be used to extrapolate the effects to larger values 
of M using the explicitly known M-dependence. Reversely, we can also use the explicit 
M-dependence to rescale all generated points to a single mass scale. We make use of such 
a rescaling in the discussion of Uckm in Section |5.4| in order to increase the statistics. 


5.1 Effective SM flavour parameters 

When dehning the ranges of the adapted flavour parameters, we have aimed at choices 
which entail a coverage of the complete range of experimentally allowed Yukawa param¬ 
eters. This turns out to be possible for most of the parameters, however, not all of them 
are distributed equally over the allowed range. 

Especially the absolute values of the CKM elements |Vtb|, |his|, iWsI and \ Vcd\ are more 
constrained in our setup. This can be traced back to fact that the direct bounds on these 
CKM elements are relatively weak while the correction AVckm (related to non-unitarity) 
is small. If we instead compare the covered range to the CKM elements as deduced from 
the SM £t assuming unitarity |^|^, we generate many points outside the 3a region. 
In particular, we End that our \Vtb\ is typically smaller. The experimental bounds we 
impose on the individual CKM elements, the ranges deduced from the SM fit and the 
ranges covered in our scan are all given in Table for comparison. Concerning the masses 
we End that they all cover the full allowed range 1^. However, there is a tendency 
to lower values especially for the charm mass. For the light quark masses, we also see a 
deEcit in the largest values of the allowed ranges, although not as pronounced as for the 
charm mass. We interpret these efiects as relics of our scan, as the setup of the model 
generally lowers the masses. The allowed ranges for the masses can be found in Table 
Having checked that the model is capable of generating the Eavour structure of the 
SM, we now turn our attention to quark Eavour efiects beyond the Standard Model. 

^In the one generation case it is possible to show that A > v^, and larger values are disfavoured in a 
perturbative theory. 
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experiment 

unitarity 

covered 


\Vud\ 

0.9736 - 0.9749 
0.9739 - 0.9747 
0.9737 - 0.9748 


iKsI 

0.2229 - 0.2277 
0.2235 - 0.2272 
0.2231 - 0.2277 


\Vub\ 

0.0027- 0.0056 
0.0031 - 0.0040 
0.0027- 0.0056 


experiment 

unitarity 

covered 


0.2010 - 0.2490 
0.2234 - 0.2271 
0.2229 - 0.2276 


IK.I 

0.9380 - 1.0340 
0.9730 - 0.9739 
0.9727 - 0.9741 


0.0372 - 0.0450 
0.0378 - 0.0450 
0.0372 - 0.0450 


experiment 

unitarity 

covered 


\Vtd\ 

0.0066 - 0.0102 
0.0079 - 0.0099 
0.0066 - 0.0102 


0.0319 - 0.0481 
0.0372 - 0.0438 
0.0355 - 0.0447 


\Vtb\ 

0.9250 - 1.1170 
0.9990 - 0.9993 
0.9711 - 0.9992 


Table 3: Coverage of the allowed range of 


(^CKM)i- 1- The hrst line corresponds to the 


direct experimental limits we impose for the scan . The second line shows the allowed 
range deduced from the SM fit assuming unitarity of the CKM matrix 
corresponds to the range covered in our scan. 


66 . The last line 


ruu 0.5-2.9MeV ma 1.2-4.8MeV 
rric 0.53 — 0.71 GeV mg 30 — 78 MeV 

mt 162 - 180 GeV mb 2.78 - 2.96 GeV 


Table 4: Ranges of the quark masses covered in the scan. 


5.2 Masses of the new heavy quarks 

The masses of the heavy up-type quark partners can be read off from the approximately 


diagonalized mass matrix of Eq. (3.7). A similar expression can be deduced for the down- 


type quark partners by replacing the index u by d. As we have explicitly factored out the 
mass scale M it is reasonable to consider only the mass ratio 


fip — 


mp 

W 


(5.1) 


With the matrix / being independent of isospin, we can directly infer that one set of heavy 
up-type partners equals one set of down-type partners in mass, denoted by U” and D’l 
in the following, where Ui {Di) are the three generations of up-type (down-type) quarks. 
This observation is numerically verihed in our scan, where the corresponding masses are 
equal at the sub-percent level. Introducing a VEV (T) ^ 0 for the SU{2) triplet flavon T 
would break this mass degeneracy by an amount proportional to (T). Being proportional 
to eu,d, the masses of the other set of heavy quark partners, denoted by 11^ and 77', will 
generally depend on isospin. 
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tan 13 


tan /S 


(a) First solution. 


(b) Second solution. 


Figure 3: Mass ratios fip of the lightest generation of heavy quark partners with respect 
to tan (3. The plots show the masses of the hrst top partner in black, the hrst bottom 
partner in blue, and the second top and bottom partners in green. The two solutions 
of Eq. (2.12) are given for each mass ratio as dashed lines in red (a) and orange (b), 
respectively. 


Focusing on the lightest generation of heavy quark partners, i.e. f, b', t" and h” 


the multiplicity of solutions to Eq. (2.11) is dominated by the multiplicity of the one 


generation case, see Eq. (2.12). The scan therefore generates two distinct bands for each 


of the four mass ratios These in turn depend on tan/3, which can be traced back 
to the appearance of yu and yd in Eq. (2.12). The resulting tan/3-dependence of the 


heavy quark masses calculated from the one generation case hts the results of the scan 
reasonably well as illustrated in Figure Here, we have plotted the mass ratios yp of the 
four lightest quark partners (t', 6', t" and b") against tan/3. For clarity we have separated 
the two possible solutions. The mass ratios of t” and b" are shown in green (as both 
are equal), the one of t' in black and that of b' in blue. Additionally, we have plotted 


the tan/3-dependence of the mass ratios deduced from Eq. (2.12) where we have used the 


central values of the SM Yukawa couplings and A = 2 as fixed input. Here, we introduce a 
colour coding for this and the following plots: masses belonging theoretically to the first 


second) solution of Eq. (2.12) are shown in red (orange) 


We can conclude from the plots of Figure that the mass ratios yp oi the heavy quark 
partners may be of order one, and thus their masses can be in the TeV-regime. However, 
depending on the solution, we may expect either a set of three quark partners (6', t" and b") 
within a small mass range or solely the top partner t' to appear in experiments. This 
shows that the two solutions differ qualitatively in their phenomenological predictions. 
We will encounter further differences in the following discussions. 


5.3 Non-unitarity of 


Due to the e\ ^-dependent contribution AVckm to the effective CKM matrix, see Eq. (3.15), 
the coupling of the bosons to the three generations of light quarks is no longer de¬ 
scribed by a unitary matrix. This non-unitarity can be quantihed by studying the uni- 
tarity triangles constructed from V^km- clarity, we limit our discussion 
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(a) The origin of the standard unitarity triangle. 
The blue-shaded regions correspond to the Icr 
and 3(7 uncertainty of the measured SM Wolfen- 


stein parameters p and rj 66 



(b) The absolute value oip + irj plotted against 
tan/3 for points with M G [1,1-2] TeV. The 
dashed blue line indicates the corresponding 


la error of \p + i'q\ as obtained in the SM 66 


Figure 4: The non-unitarity of V^km- For the hrst (second) solution of Eq. (2.12), we 
show p + irj, i.e. the measure of the deviation from unitarity, in red (orange). 


to the “standard unitarity triangle”. However, we have checked that the other unitarity 
triangles lead to analogous results with similar or smaller effects. 

Generically we do not expect large deviations from the SM as we have limited our scan 
to points which are compatible with the absolute values of the CKM matrix elements, 
taking into account the uncertainty of the measurements. This is a simplifying assumption 
as the determination of the CKM elements from experimental measurement is also affected 
by the couplings of the right-chiral quarks. To give a complete picture, we would have 
to redo the extraction of the CKM elements in the presence of additional couplings of 
right-chiral quarks (Cckm) and without the assumption of unitarity. Clearly, such an 
explicit £t is beyond the scope of our current work. 

In order to quantify the non-unitarity of via the standard unitarity triangle, we 
dehne the complex quantity 


p + ip = 


VubV:, + v,kV:^ + Vt,v, 






Vr,V* 


( 6 . 2 ) 


whose deviation from zero is a measure of the deviation from unitarity. In Figure]^ we 
show 'p + irj for the complete set of viable points at the origin of the unitarity triangle. In 
addition, we present the tan/3-dependence of its absolute value [p + irj\, where we limit 
ourselves to viable points which lie in the mass band M G [1,1-2] TeV for the sake of clarity. 
The points associated with the two solutions are again distinguished by the colour code 
introduced above. We compare these results of our scan with the la and 3 ct uncertainties 
of the measured SM Wolfenstein parameters p and p 66 . As expected, the deviation from 


unitarity predicted in our scan is within the SM uncertainties. Moreover, we again find 
that the size of this effect depends crucially on the type of solution. This is particularly 
apparent in the tan/3-dependence of \p + ir]\, cf. Figure 4(b) Here, the second (orange) 


solution is basically independent of tan /3 whereas the hrst (red) is of similar size for 
tan/3 = 1, but becomes negligible for tan/3 > 3. Thus, the effect can be practically 
absent even for small values of the Favour breaking mass scale M. Considering the values 
oi p + irj in the complex plane near the origin of the unitarity triangle, we observe that 
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Re(Uub/ 


(a) Impact of Uub on the determination of \Vub\ 
as given in The coloured bands show dif¬ 
ferent measurements: inclusive decays (blue), 
B — )■ ttIv (red), B — >■ plv (yellow), B ^ tv 
(green). The points obtained in our scan are 
shown as black dots. 



Re(Ucb/Vcb) 


(b) Impact of Ucb on the determination of 11461 
as given in [^. The coloured bands show dif¬ 
ferent measurements: inclusive decays (blue), 
B —)■ D*lv (red), B —>■ Dlv (yellow). The points 
obtained in our scan are shown as black dots. 


Figure 5: The impact of Uqkm on the determination of |Vckm|- Using the explicit M- 
dependence, the points of our scan have been scaled to the reference scale M = 1 TeV. 


the real part p is generically positive, cf. Figure 4(a) This indicates that the effect is 


dominated by the decrease of \ Vtb\ which in turn reduces the length of the upper right line 
of the unitarity triangle. 


5.4 Coupling of right-chiral quarks to via Uckm 


As described in Section |3.2[ the light generations of right-chiral quarks couple to the 
SM fU^ gauge boson due to the mixing of qr with Er in Eq. (2.4). The corresponding 


coupling strength is parameterized by Uckm as stated explicitly in Eq. (3.15c). Formally 
it is proportional to which we hnd conhrmed in our scan. Furthermore, we see no 


direct correlation between Uckm and Vckm- The existence of such a coupling modihes the 
determination of the CKM elements from experimental data. In particular, it might be 
signihcant in view of the puzzling tension between the determination of \ Vub\ and \ Vcb\ from 
exclusive and inclusive U-meson decays within the SM 66,70 . The additional couplings 


Uub and Ucb of right-chiral quarks to can in principle reduce (though not completely 
resolve) this tension as the inclusive decay is proportional to |14bp -b lUxbl"^, while the 
exclusive decay is proportional tcgii4b±u,b|2, see e.g. the discussions in (7l}j7^ . Yet, from 
the point of view of our model, our scan may reveal which experimental measurements 
are “theoretically preferred”. 

In our scan, we have searched for such effects, but do not hnd a conclusive answer. 
Figureshows the la regions of \ Vub\ and \ Vcb\ as a function of Uub and Ucb 72 . We have 


additionally plotted the 3,000 points of our scan over the full parameter range, scaled to 
the reference mass scale M = 1 TeV. As a result we hnd that the relative contributions 
of Uckm are at the percent level at most. Furthermore, the scan shows no preference for 
any of the measurements, and no apparent correlation between Uub and Ucb is seen. 


^Here the sign is determined by the currents (axial or vector) mediating the decay. 
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(a) Correlation between the coupling of the Z 
boson to the left- and right-chiral bottom quark. 
The dashed grey line indicates the diagonal. 



(b) The dependence of on tan/3. The 

dashed blue line indicates the Scr limit of Eq. (5.3) 
(points below are allowed). 


Figure 6 : Correction of the Z boson coupling to the b quark. In red (orange) we show 


the points corresponding to the hrst (second) solution of Eq. (2.12) for the mass band 
M e [1,1.2] TeV. 


5.5 Anomalous Z coupling 


As can be seen explicitly in Eq. (3.14), the coupling of the Z boson to the left- and right- 


chiral quarks is modihed as well. Similarly to the effects discussed above, the largest 
deviation arises in the top quark coupling, where it is in the range of a few percent for 
moderate values of M. Experimentally, the Z-top-top coupling cannot be determined 
directly from the decay of an on-shell Z boson as it is kinematically forbidden. Thus 
top physics processes such as ttZ production at the LHC 74 , or the production of a 
pair of top quarks mediated by the neutral gauge bosons Z* and 7 * at the ILC 75 
provide alternative and promising ways of measuring the coupling of the Z to the top 
quark. However, as the resulting bounds are expected to be relatively weak, we focus on 
the coupling of the Z to the bottom quark. This decay has been measured precisely at 


LEP2 76 with r(Z —>■ bb) ~ 375.87± 0.17 MeV 66 . It involves the couplings of both the 
left- as well as the right-chiral b qua rks. C onsidering the correlation of their corrections 

we observe that 


within our model, plotted in Figure 6 (a 


and 


^9zbnbn\ 


roughly of equal size, unless the values are small and thus negligible for phenomenological 
purposes. We can therefore make the simplifying assumption that both corrections are 
approximately identical. Assuming furthermore that the correction to the Zbb coupling is 
within the experimental uncertainty, we can set the following 3a limit on the couplings. 


\^9zb 


ZbLbL I 


\^9, 


ZbRbR 


< 6.8 X 10 


-4 


(5.3) 


Fig. 6 (b) 
(orange) 


This bound reduces the number of viable points of the scan by roughly 10 percent. The 
so-e xclud ed points all belong to the hrst solution of Eq. (2.12). Assuming M G [1,1.2] TeV, 
shows the tan/3-dependence of for both solutions. While the second 

solution does not vary with tan/3, the hrst (red) is characterized by decreasing 
values of [A/z^j^^^l for increasing tan/3. Thus, the bound of Eq. (5.3) is satished more 


often for larger values of tan/3, and all model points corresponding to the hrst solution 
with tan /3 < 4 can be excluded (yet smaller values are allowed for larger M). 
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(a) The normalized coupling as a function 
of M. The SM expectation is shown by the 
dashed blue line. The error of the experimental 
direct measurement exceeds the plotted range. 



I^ibl 


(b) Correlation of the normalized coupling 
and \Vtb\- The point corresponding to the SM 
value is shown in blue. The experimental errors 
on both quantities exceed the plotted ranges. 


Figure 7: Correction of the normalized Higgs-top-top coupling In red (orange) we 
show the points corresponding to the hrst (second) solution of Eq. (2.12). 


5.6 Anomalous Higgs coupling 


Following the discussion of Section |3.3[ the Higgs coupling to the top quark can 
receive sizeable corrections in our setup. However, such an effect is negligible already 


for the b quark as well as for all lighter quarks. Being of order e^, see Eq. (|3.17), the 


correction to reduces drastically for increasing M. In order to visualize the deviation 
from the SM, it is useful to normalize the coupling g/^j^ to its SM expectation \/2mt/vu- 
To this end, we define 

Ihtt = - 7 ^ - • (5.4) 

V 2 mt 

shows that the deviation of this quantity from its SM value is 


The plot in Figure 7(a 


of the order of a few percent for small values of M. Experimentally, the Higgs-top-top 
coupling is however only poorly known with an uncertainty of about 30 percent 77 79 


The correction to is therefore expected to be more important for precision observables, 
where the top quark and the Higgs appear in loops so that indirect constraints apply. 
Searching for correlations, we find that and \ Vtb\ show a similar M-dependence for 


the points corresponding to the second (orange) solution of Eq. (2.12), cf. Figure 7(b) 


Such a correlation can be of phenomenological interest as the top quark decays nearly 
instantaneously into a bottom quark and a W boson. For smaller values of M, it thus 
amplifies the reduction of associated Higgs production, gg — ttH, a process currently 
searched for at the LHC 80,81 . In contrast to this, the first (red) solution of Eq. (2.12) 


does not feature such a correlation. 

Besides the reduction of the normalized coupling 'jhqq, we also obtain flavour changing 
couplings of the Higgs boson which can be relevant for flavour changing neutral currents 
(FCNC). The corresponding weak experimental bounds for heavy quarks will eventually 
become stricter. A detailed analysis of FCNC effects, including the Z boson coupling, is 
therefore well motivated but beyond the scope of this paper. Concerning the light quarks, 
for which strict bounds on FCNC exist, we stress that these are naturally suppressed in 
our PS flavour model. 
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6 Conclusion 


Combining the idea of vertical and horizontal unihcation, we have constructed a Pati- 
Salam symmetric model with gauged SU{3)j x SU{3)i]: flavour symmetry. In the un¬ 
broken phase, the model is characterized by an explicit left-right (i.e. Z 2 ) symmetry and 
an electroweak symmetry breaking Higgs bi-doublet. Analogous to the construction of 
Grinstein, Redi and Villadoro (GRV), the flavour structure of the SM quarks and leptons 
originates in a renormalizable Lagrangian involving new heavy fermions as well as a set 
of scalar flavon helds. Vacuum expectation values of the latter generate masses for the 
additional fermions, which in turn mediate the associated flavour symmetry breaking to 
the SM fermions in a seesaw-like fashion. In contrast to GRV, our left-right symmetric 
PS setup requires the introduction of two (rather than one) heavy fermionic partners for 
each SM fermion. As a consequence, the cancellation of flavour gauge anomalies entails a 
further extension of the fermionic particle content. This hnds a natural realization in the 
lepton sector where light Majorana neutrino masses can be generated by introducing PS 
neutral fermions which acquire large Majorana masses from another set of flavon helds. 
These havons break the havour symmetry completely at a very high scale so that all 
havour gauge bosons decouple from low-energy physics. 

In the quark sector, the effective SM Yukawa matrices are related to the VEVs of the 
havon helds neither linearly as in the MFV approach nor in a simple inverse way as in the 
original GRV setup. The approximate relations, obtained by integrating out the heavy 
fermions and given in Eq. (2.11), do not hx the havon VEVs uniquely. Besides the un¬ 
known mixing of the right-chiral quarks (which is physical in our PS model), there exists a 
total of eight solutions for the havon VEVs, all giving the same approximate SM Yukawa 
matrices. Inserting these vacua into the complete fermion mass matrices, it is possible to 
extract the masses and the mixing angles of the three generations of SM quarks. Gompar- 
ing these with the physical values, we hnd that the approximate formulas of Eq. (2.11) do 
not adequately describe the third generation. Such a behaviour is expected and can be 
traced back to substantial mixing of the third generation of quarks with the corresponding 
heavy fermionic partners. In order to correct for this ehect, we dehne adapted havour 
parameters and take these as input for our numerical scan. Similarly to the GRV case, 
the construction of the model guarantees that non-standard havour transitions between 
the hrst and the second generation are highly suppressed, while new physics ehects in¬ 
volving the third generation of quarks are expected to be phenomenologically relevant. 
Gharacteristic ehects include the non-unitarity of the ehective GKM matrix, the coupling 
of the SM right-chiral quarks to the electroweak boson and anomalous couplings of 
the Z as well as the Higgs boson to two SM quarks. Quantitatively, all these ehects are of 
order where eu^ = {h^^)/M denotes the ratio of the Higgs VEV over the sole explicit 
mass scale M of the model. By construction, M corresponds to the scale of the lightest 
fermionic partners, and we choose it to be in the TeV regime. The discussed new physics 
signatures of our PS GRV setup are therefore below the current experimental sensitiv¬ 
ity. Yet, depending on the type of solution for the havon VEVs, the phenomenological 
predictions beyond the SM might be tested in the near future. 

In this work we have not addressed the question of the underlying dynamics which 
governs the vacuum structure of the havon helds. In our numerical scan we hnd that 
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viable flavon VEVs feature hierarchical patterns which are more or less comparable to 
those of the inverse Yukawa matrices. Ideas of generating hierarchical flavon VEVs from an 
appropriate flavour-invariant scalar potential have already been discussed in the context of 
MFV, and we intend to apply similar methods of spontaneous flavour symmetry breaking 
to our PS GRV setup in a future publication. Another subject of further investigation 
concerns the lepton sector. In this work, we have presented a neutrino extension in 
which the flavour structure of the neutrinos decouples from the hierarchical patterns of 
the quark sector. The construction of a fully realistic charged lepton sector, including 
a detailed discussion of lepton-flavour violating processes and lepton non-universality, is 
left for future work. 
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Appendix 


A Lepton sector 

A.l Neutrinos 


In order to accommodate Majorana neutrinos, we have extended the Pati-Salam GRV 
model by /lyuk of Eq. (2.5). Integrating out the Majorana fermions 0^ and Qr generates 
an effective heavy Majorana mass for cf. Eq. (2.13). This induces a Majorana mass 
term for the left-chiral neutrino ql via the terms of Eq. (2.4). In the hrst part of this 
appendix, we derive the structure of the resulting neutrino mass matrix. In the second 
part, we discuss an alternative ansatz in which (rather than acquires a heavy 
Majorana mass. This is one of several alternative examples leading to light neutrino 
masss via Eq. (2.4). The havour structure of the resulting mass matrix is, however, less 
attractive as will be discussed exemplarily below. 


A. 1.1 Preferred neutrino extension 


Starting with the neutral part of the Lagrangian of Eq. (2.3), it is straightforward to 
integrate out the heavy helds Ql,r and El.r- Inserting the flavon and Higgs VEVs we 
hnd 


Cass = + + + h.c. 

-A^gls-^gj^ +h.c. 

- +h.c., (A.l) 
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with Myv = 


_ 0-1 


cf. Eq. (2.13). To arrive at this simple intermediate result, we 


have assumed ^ Next, we integrate out to obtain 


O' 




2 + M j s.[X^^ql + M j + h.c. 

+ M - A ^ 5"^ g^ + h.c., (A.2) 

where we have separated the Majorana mass term in the first line from the Dirac terms in 
the second line. and g^ form a Dirac pair and can be integrated out simultaneously. 
It is interesting to note that already the equation of motion of g^ is sufficient to deduce 
the effective light neutrino mass term. Inserting —)■ A g^ yields 

2{ipCKY (l + s 


L 


V 

mass 


— 2^L 




(A.3) 


which corresponds to the effective mass matrix reported in Eq. ( |2.14 ). Following this 
derivation, the occurrence of the second term in the sum (1 + s~^tu) is traced back to the 


coupling of to g^ via E^, and g^, cf. Eq. (A.l). 

It is important to note that a high scale VEV {SD is essential. A vanishing or too 
small value would imply only a small Majorana mass for 0^ as well as a mixing with E^. 
Integrating out the heavy fields E^, E^ and g^ results in a strongly hierarchical Dirac 
mass term for g^ and 0^ which is phenomenologically excluded, unless the Majorana 


contribution, identical to Eq. (2.14), to the neutrino mass dominates. 


A.1.2 An alternative neutrino extension 

Altering the flavour quantum numbers of $ & from (8,1) & (1,8) to (3,3) & (3,3), 


entails a slightly modified form of /lyuk Eq. (2.5): the Z 2 pair simply gets 

replaced by {q^, qn). Analogous to the preferred case, one obtains a Majorana mass term 
for g^ after integrating out Ql,r- 


Cass = Ag^^E^ + ME^f„E^ + ME^g^ + h.c. 

+h.c. 

~ 2 ^R ^'r + ’ 


(A.4) 




where ip''^ p'. The heavy right-chiral neutrino g^ couples to g^ both directly 

and indirectly via E^ and E^, thereby generating an effective Majorana mass for the left- 


chiral neutrino g^. Integrating out g^ and E^^ in Eq. (A.4) results in 


C 


u 

mass 


= -n^ 
2 




1-1 


{K^ As ^) 


2A2 






(A.5) 


Comparing Eq. (A.5) and Eq. (A.3) we notice two important differences. First, due to 


its transformation as a (3,3) in flavour space, the matrix p' contributes non-trivially to 
the flavour structure of the light neutrinos. Second, the factor {t~^ + s 




Yu features 

strong hierarchies in contrast to (1 -|- s~^tu). These hierarchies must be approximately 
compensated by an even stronger hierarchy in in order to generate a phenomenolog¬ 
ically viable neutrino mass spectrum. As such a cancellation of hierarchies appears rather 
ad hoc, we abandon this alternative neutrino extension. 
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A.2 Charged leptons 


As already mentioned in Section |2.3| , the structure of the charged lepton Yukawa matrix 
is identical to the one of the down-type quarks. At the GUT scale we have Y^. = which 
is a reasonable approximation. Deviations from this equality can be induced by allowing 
the flavon helds S or T^''^ to transform non-trivially under SU{A). Due to their appearance 


in the Lagrangian of Eq. (2.4), the only alternative representation is the adjoint 15 of 


SU{A). Such a choice does not alter the quark sector, but it introduces a Georgi-Jarlskog 
factor of —3 for the leptons [^. For instance, if S is kept in the singlet representation of 
SU{A) while T^''^ furnishes the adjoint representation, we obtain a difference between the 
Yukawa matrices of the charged leptons and down-type quarks. Explicitly, and adopting 
the approximation of Section m we obtain 


Y, 


{s -1') 


^ + s-^ 




(s -|- 3 


+ s-^ 


(A.6) 


More complicated relations are possible if we allow the flavons S and/or to appear in 
both SU{Y) representations. We conclude that the model can therefore accommodate a 
viable charged lepton sector as well. 

B Determining s and t' numerically 


As already discussed in Section 2^, it is a not-trivial task to invert the relations between 
the flavour breaking VEVs and the Yukawa matrices. In this Appendix, we present a 
procedure for calculating s and t' numerically for a given pair of and Y^ assuming 


Eq. (2.11) to hold exactly. We therefore write 


— jYu — {s 1') ^ -)- s 


Dehning H = s we obtain the expressions 
-}-YuS = il + H)-^ + l, 


-lY,= {s-t')-^ + s-^ 


lYdS = {l-H)-^ + l 


which can be combined to remove the s-dependence 

G = r,y-‘ = [(1 - +1] [(1 + i/)-‘+1] 


-1 


(B.l) 


(B.2) 


(B.3) 


Let us assume for the moment that G is a triangular matrix, which can always be achieved 
by a suitable similarity transformation. In that case, also H has a triangular form, so that 


Eq. (B.3) can be solved explicitly for the six elements of H. We thus obtain an explicit 


solution for H{G) in the special case of a triangular matrix G. 

As Yd Y~^ is generally not given in a triangular form, we need to generalize the solution 
for H (G) to generic forms of G. For this purpose we expand the function H (G) as a series 
in the matrix G. Due to the Gayley-Hamilton theorem this series stops after the quadratic 
termj^ Hence, we can express H{G) in a basis independent way by 

H{G) = al + bG + cG‘^ , (B.4) 


®In general the theorem states that for each n x n matrix M, M" can be expressed as a polynomial 
J2i<n where the factors Xi depend solely on traces of powers of M up to order M”. For a detailed 


discussion see e.g. 82 83 . 
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where the coefficients a, b, and c depend only on the traces of G, and G^. In the 


special case of a triangnlar matrix G, it is possible to solve Eq. (B.4) analytically for a, b 


and c. As these coefficients are basis independent, Eq. (B.4) holds trne in any basis; in 
particnlar G = Y~^ need not be of triangnlar form. 

In principle, we can determine H for a given pair of Yukawa matrices Y^^d analyti¬ 
cally. Due to their enormous length, the resulting expressions are however not particu¬ 
larly instructive which is why we content ourselves with a numerical evaluation. Having 


determined the matrix H, we can use Eq. (B.2) to calculate s 


s = -AF„-‘((! + //)-'+ 1) , 


(B.6) 


and the dehnition of H to determine t' = sH. Thus, we have obtained numerical matrices 
s and t' for a given pair of Yukawa matrices Yu^d- 


C Diagonalizing the quark mass matrices 

In this appendix, we present the technical details related to the diagonalization of the 


quark mass matrices, discussed in Section 3.1 Starting from the original basis given 
in Eqs. ( [T2p3| ), we perform a sequence of basis transformations to diagonalize JYY to 
second order in e„. The individual intermediate bases are labelled by a subscript (z). 

1. Basis with diagonal s = s: 

Thanks to the flavour symmetry SU{3)i x SU{3)n we can always choose a basis in 


which the matrix s is diagonal. is therefore identical to of Eq. (3.3) with 

s —)■ s. Note that this choice of basis is valid for both the up and the down sector. 

Diagonalizing 

In the second step, we apply the basis transformation 

¥1,2, = Till, diag(l,v;t,l) , n( 2 ) = diag(K.,(7.,l) n(i) . (C.l) 


to render diagonal, see Eq. (3.5). The resulting mass matrix takes the form 

UlXe„, 


-M(2) = 



i 


M . 


(C.2) 


Notice that the unitary rotations Vu and Uu reappear at order only. It is further¬ 
more important to realize that the corresponding matrices Vd and Ud for the down 
sector differ from those of the up sector. 

3. Diagonalizing Af“ 2 ) for Cu = 0: 

Setting e„ to zero, the mixing of the three generations disappears completely. In¬ 
troducing the following cosines and sines. 


tl 




\/l + 


c, = 


l + tl^ 


s, = 


, (C.3) 


l + tl^ 
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we can define the diagonal 3x3 matrices 


Cx = diag(4,4,4) , = diag(s^,4,s^) , (C.4) 

where the index x stands for s oitu. With this notation at hand, it is straightforward 
to diagonalize M.^ 2 ) ~ basis of 


-^u -^u 


/ Cs 

0 

■sA 

1 j 

fct^ 

^tu 


0 

1 

0 

1 ’ ^^(3) - 1 

^tu 


0 

V-Ss 

0 

cj 

^0 

0 

1/ 


R(2) 


(C.5) 


the fnll mass matrix, inclnding the terms proportional to e„, is given by 


■^(3) = 


( CsU^st^)X6u ( + Csf/^Q^)A 0 

. . ° „ S' . I ^ 

(csV^Ct„ 6u + SgU^Ct^^Xeu s 


-1 


Identifying this with the matrix of Eq. (3.6), dehnes the 3x3 matrices a^, bu, Cu, 
du, Cu, and /. The latter two are diagonal and positive dehnite. 

4. Block-diagonalizing Wl“ 3 ) up to order e^: 

Having parameterized -Ad^s) as in Eq. (3.6), we proceed to block-diagonalize this 
matrix np to order e^. To this end, we adopt the following nnitary transformations. 


= 'I’m |Ki 2 (^r 2 )|' |K 2 s(a)]' [Ki3(ff3)]* . (C.7a) 

'tfl,!) = lKl2(Cr2)l 1K23(C2“3)1 |Ki3(Ci“3)1 'I'r( 3) > (C-^b) 

where TZa/siO denotes a “rotation in the a-(3 plane”, expanded to second order in 
For {a,(3) = (1,2), one has 

-e o\ 

7^l2(0= l-l^e 0 , (C.8) 

\ 0 0 1 / 


and the expressions for the other two pairs, (2,3) and (1,3), are identical np to 
obvious permutations of rows and columns. The matrices ^ and ( are given in terms 
of the parameters of Eq. (3.6) by 


— h f: IfUl _ __ 

S12 — , LS23jp ~ 2 f 2^'^ 

~ Jj 


Ci2 = alKe. 


-2 2 
u u^u ’ 


[C23]i7 ~ ..2 




er - fj 


313 


313 


aucif ^ el , 

(C.9a) 

^uJ • 

(C.9b) 


We note that the matrices ^23 C 23 cannot be written as a simple product of 

matrices. With a little bit of algebra, it is possible to show that the resulting mass 
matrix is given by the matrix in Eq. (3.7). 
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5. The approximate mass basis: 


The hnal step of the sequence of basis transformations consists in the diagonalization 
of the upper left 3x3 block of -^“ 4 ), see Eq. (3.8). The corresponding transformation 
reads 


diag(Vi,l,l) , ni,) = ■ (C.IO) 

The basis corresponds to the approximate mass basis denoted by in 

Section |3l 


Applying this sequ ence of basis transformations to the gauge-kinetic terms changes the 
/C matrices of Eq. (3.13) to 


r+ — 

^L(5) - 


( v„- V„ i + ^Ut2) v] v„ ^^2 -V„ (era - \ 


pu tilt 
S12 


\ (?ra^-efa^+e 2 VCr 2 ^)v, 

u^CUizUl 


Cu ^ Cd I Cd ’ 
S12 S12 “t" S23S23 
— pd t 
C2a 


^23 

^ ~ 5 (^ 2 a^C 2 a + C^a^Cla) / 


+ ^(4,d) ) 

(C.lla) 


r+ — 

'^R{5) — 


iiasia '^‘d ^«si 3 S 2 a 

AU t ; yt J-U /-d t 

S 2 asia ‘^d ‘52aS23 

pd ^ 1 A Ad t 

Sia ‘^d S 23 


-Uu era 

-C 2“3 I + ^(4,d) , 

1 (pu t pu I pu t pu I pd t Ad I Ad t Ad 1 
2v‘5l3 S13 “I" S23 S23 Sl3 Sl3 S23 S237 


/vei 2 ei 2 ^vt 
- 62 ^ vt 
\—^23^^12^ V'l' 


— V ^12 

^23^ 


Al( 5) — I —^12^^^ 1 — ^12^^12 — ^23^23^ ^23 ) + d’(e^) 


— E ^12^23^ 
^23 1 ' 

^23^^23 / 


(C.llb) 

(C.llc) 


y'O _ 

^i?(5) - 


CiaCia^ ^ Ci3C23^ 

C23Cl3^^^ C23C23^ 

~Cl3^^^ ~C23^ 


Cl3 
— C 23 

1 — Cl3^Cl3 + C23^C23t 


+ 0{e^) 


(C.lld) 


As the basis transformations depend on the isospin components, we have labelled the re¬ 
sulting flavour structures by superscripts which indicate whether they relate to charged or 
neutral currents. In the latter case, we have suppressed the isospin indices for convenience. 
Note that can be derived from by dropping the indices u and d. 
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